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Q\ • Abstract 

Starting from the QCD Lagrangian and taking into account both 
Q^l perturbative and nonperturbative effects, we use the method of vac- 

o ■ uum correlators to derive the Dirac equation (rigorously for the Coulomb 

D . interaction and heuristically for the confining potential) for a system 

consisting of a light quark and heavy antiquark. As a result the con- 
fining potential is a Lorentz scalar, and the Coulomb part the fourth 
^ ■ component of a 4-vector. The energy spectrum of the Dirac equation 

H I is considered for these potentials. Numerical calculations of energy 

eigenvalues E = Enn are performed, and some exact solutions of the 
Dirac equation in the case E = are found. An effective-potential 
method convenient for qualitative study of the solutions of the Dirac 
equation is developed. The connection with experimental spectra of 
D- and 5-mesons is briefly discussed. 

1 Introduction 

The quark-antiquark systems consisting of one heavy quark (antiquark) 
Q and one light antiquark (quark) q are QCD analogs of the hydrogen atom 
and thus are of fundamental importance. 
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Recently the issue of the heavy-hght bosons has become a topic of vivid 
interest in both analytic and Monte Carlo QCD studies [1]. 

From the theoretical viewpoint the interest in heavy-light systems stems 
from several considerations. First, in the limit of one infinitely heavy quark, 
one hopes to get the dynamics of a light quark in the external field of a heavy 
one. That would be similar to the picture of the hydrogen atom. 

Second, since the external field is time- independent, one may hope to 
obtain a static potential in QCD together with spin-dependent forces, as has 
been done for heavy quarkonia [2]. 

An important issue in this connection is the Lorentz nature of the con- 
fining part of the potential. Arguments in favor of a scalar nature are found 
in the form of inequalities [3] and also in the form of spin-dependent terms 

[2,4]. 

Third, in the heavy-light system one may study how the chiral symmetry 
breaking (CSB) affects the spectrum. When one quark mass is vanishing, in 
the chiral symmetric case the spectrum would consist of parity doublets, and 
the CSB would hft the degeneracy. 

Fourth, using the Dirac equation we implement explicitly relativistic dy- 
namics and can study relativistic properties of the spectrum, e.g., in the 
case of a vanishing quark mass, and also relativistic spin splittings in the 
spectrum. 

In particular, in our case the spin symmetry now being widely discussed 
[1] is present, since the spin of the heavy quark is decoupled. Hence every 
state of the Dirac equation with a given j and parity corresponds to two 
almost degenerate states of the (qQ) system with j and the heavy quark spin 
S = 1/2 adding to J = j ± 1/2. 

Our final point is that the Dirac equation description yields a dynam- 
ical framework for the heavy-light mesons which can be used to calculate 
meson matrix elements and form factors to compare with experiment and 
semiphenomenological approaches now widely used in this context [1,5]. 

Here we attempt to derive explicitly the Dirac equation for the heavy- 
light quark system starting from the QCD Lagrangian and incorporating both 
perturbative and nonperturbative effects in the framework of the vacuum 
correlator method (VCM) [6-8]. 

In the course of the derivation we use some approximations, which are 
clearly stated at each step, e.g., neglect of virtual quark pairs (the quench- 
ing approximation), keeping the lowest (quadratic) cumulant in the cluster 



expansion, and the "limit of local dynamics" [2], Tg — ^ 0, where Tg is the 
gluonic vacuum correlation length. We discuss physical implications of these 
approximations in the main text below. 

Even with those approximations we are unable to prove rigorously the ap- 
pearance of a static scalar confining potential, but we give strong arguments 
in favor of it. The situation is much better for the perturbative part, and 
the existence of a external vector potential in the limit of one heavy quark 
mass is explicitly shown. 

At this point we start with the Dirac equation containing a vector Coulomb 
part and a scalar confining part, and we study properties of its spectrum. 
For comparison we also consider two other cases: i)the equation has a vector 
confining part and a Coulomb part; ii) both parts of the potential transform 
as scalars. We are explicitly interested in the limit of a vanishing light quark 
mass and the spin dependence of the energy eigenvalues. 

We show explicitly that a reasonable spectrum occurs only for a scalar 
confining part. In this case the scalar potential breaks explicitly chiral sym- 
metry, and states with opposite parities are not degenerate. In the case of a 
vector confining part, only quasistationary states are found. 

We write possible general symmetry properties of spectral states and also 
the corresponding properties in the particular cases of a zero scalar potential, 
a zero vector potential, and a zero mass. 

We also find some exact analytic solutions of the Dirac equation which 
occur for i? = 0. This enables us to find an equation, which determines 
that value of the Coulomb constant {( = (nx), for which the (njl) discrete 
energy level reaches the value E = 0. We note an analogy of this problem 
with the relativistic Coulomb problem of an electron in the field of a heavy 
supercritical nucleus with Z > 137 [9-11]. 

The paper is organized as follow. In Sec. 2 we use the Feynman- 
Schwinger representation of the quark-antiquark Green function to derive 
the limit of one infinit mass. 

In Sec. 3 we compare this limit with the Dirac equation and discuss the 
Lorentz nature of the confining interaction. 

In Sec. 4 we discuss the properties of the spectrum of the Dirac equation. 

In Sec. 5 we discuss symmetry properties of the spectrum, in particular 
chiral properties. 

In Sec. 6 we discuss numerical results for the spectrum of the Dirac 
equation in the limit of a vanishing small quark mass. 



In Sec. 7 we investigate the exact analytic solutions of the Dirac equation 
which occur for E = 0. 

In Sec. 8 we summarize our results and compare the Dirac spectrum 
obtained here with experiment and other approaches. This enables us to 
also discuss large-mass corrections. 

2 Feynman— Schwinger representation for the 
heavy— light Green's function 

Consider the quark-antiquark system of one heavy antiquark (with mass 
7712) and a light quark of much smaller mass rrii, which we neglect in some 
cases below. 

The Green's function G of such a system, with the quarks initially at 
points y, y and finally at points x, x, is given by the path integral over 
fermionic and gluonic fields A in ref. [8]. Evaluating the fermionic integrals 
and neglecting the fermionic determinants (the quenching approximation) for 
simplicity, one arrives at the amplitude depicted graphically in Fig.l, where 
Aj is the quark Green's function with quark mass rrii, and the F's belong to 
the initial and final wave functions: 

with a similar expression for \I'(x, x). Thus we can write 

G{xx\yy) = (tr[r(x, a;)Ai(a;, y)T{y, y)A2{y, x)]). (1) 

The angular brackets here denote the integration over gluonic fields. This 
integration we perform using the cluster expansion [6,7]. We arrive at the 
following expression involving proper-time and path integrals for the quark 
{s,z) and the antiquark {s,z): 



00 00 
G{xx\yy) = ds dsDzDze 



-K-K^ 





xtr[r,(mi - ^(i))7,,-rj,(m2 - D^^% (2) 



where 



1 f 
K = mis + - z'^dr, 



K = m\s H — j z df. 



In (2), r^, Ty refer to the Lorentz structures in the initial and final states, and 
D^*) and 'jzz are given in ref. [7]. We rewrite the latter here in the following 
form (in the lowest order of a cluster expansion, where perturbative and 
nonperturbative contributions factorize [8]): 



725 = exp 
where 



g' 



df^v{u)dfpx{u){F^y{u)Fpx{u')) 



(3) 



df^,u{u) = ds^^{u) - i(J^^}dT + iaf^dr, (4) 

and {F F) is the gluonic correlator. The latter can be split into perturbative 
and nonperturbative parts as in ref. [7]. 

Now we discuss the limit of one heavy mass, m2 -^ oo. In this case, when 
1712 is much larger than the interaction strength, particle 2 is moving along 
the straight-line trajectory 



Zi, = ypA = r, < r < s, 

s 

XA-yA = T, — = 2m2, dr = ; (5) 

5 ZTfl2 

Xi = Vi = Ri is the position in space of particle 2. 

One can see that in this limit the spin interaction of particle 2 (terms 
a'^'^^df) is ~ 0(1/1712) and can be neglected. 

The perturbative part of the cumulant {{F{x)F{y))) can be written as 

[6,7] 

g"^ / da^^{u) / da px{u){{Fpy{u)Fpx{u)))pert = 



c c 



g'^ / dzp / dz' {Ap{z)Ap{z'))pert 



C2g' f , f dz'^ 

(z — z 



4^,2 / ^^^ I f._l,2 ^ (6) 



c c 



where C2 is the quadratic Casimir operator, C2{Nc = 3) = 4/3, and C is 
the closed contour depicted in Fig.l. The integral in (6) is divergent and 
is regularized by inserting a Z factor in front of the exponential in (3) and 
introducing the minimal distance 6 of the points z and z' in (6) (for details 
see [12]). 

Therefore we can now study the situation when only the perturbative 
interaction is present, and the Green's function looks like 

00 / '^ \ 

G{xx\yy) = / dsDzexp —mis — - z^dr tr[ra;(mi H — z) x 



xZexp 



/ / ids^„{u) + -(jl^^dT)dspxiu){{Ff,^{u)Fpxiu')))pert 



xr,(m2(l-74))], (7) 



where we have neglected terms representing the self-interaction of particle 1. 
Our goal now is to rewrite the exponential function in (7) in the form 



exp 



/■ - 1 - 

W J {A^{z)dZf, + -a^^^F^ydr) 



(8) 



(a similar derivation for QED was done in [13]), where we have defined (fol- 
lowing [6]) 



A^{z) = ig I {A^{z)Ap{z'))dz'p 



47r I z - R 



= ^ 

47r2 

■C2. 



dz' 



[z — Z' 



'\2 



Similarly, 



flU 



u^Ai, — dyAp. 



(9) 



(10) 



One can also verify that the explicit expression for {{F F))pert, which one 
obtains from the Feynman gauge propagator [6], leads to 



{Ap{z)A^{z'))pert 



5, 



fiu 



47T^(Z-Z'V' 



(11) 



Hence we can write the qq Green's function, keeping only the perturbative 
interaction: 

G{xx\yy) = tr[r^G'(a;, y)Tym2il - 74)], (12) 
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where we have defined 

G{x,y) = I dsDz[mi - D{A)]x 

(S X s , 

-mjs -J z^dr + ig A^dz^ + g a^^F^^dr 
y ' 

f ds[mi - D{A)]exp {-s[m^ - D^{A)]} \ y) 



y 

OD 

[X ' 




= {x\ {m, + D{A))-'\y). (13) 

We observe that G{x, y) satisfies the Dirac equation 

K + b{A)]G{x, y) = d^'^Xx - y). (14) 

In the next section we study the effect of the nonperturbative interaction. 

3 Confining force and the Dirac equation 

We now turn to the nonperturbative interaction, which provides confine- 
ment due to the presence of the special (Kronecker) structure in the quadratic 
cumulant [6]. Again neglecting the self-interaction of light particle, we find 
an equation of the same form as (7), but we should add {{F F))pert also the 
nonperturbative part {{FF))nonpert, which can be written as [6] 

{{F^,{u)F,x{u'))) = {6^,6,x - S^xSup)D{u - u') + 

+ 2 i Q — [(« - u')p5^x -{u- u)x5^,p] + ^p^v\\ Di{u - u'). (15) 

It has been shown [6,7] that only D (and not Di) yields confinement (the 
area law of the Wilson loop). We concentrate first on this term, disregarding 
also the spin term a^^^} in (7). We have 

ds^^{u) J dspx{u'){{Fp^{u)Fpx{u'))) -^ J ds^^{u)ds^^{u')D{u - u') + ... . 

(16) 



Using the same arguments as in [6], one can see that we get the area-law 
exponential function in (7). Namely, for large area we have 

exp{-aSmin), ^ = 7; d^uD{u), (17) 
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where Smm is the minimal area inside the contour formed by the straight-line 
trajectory of the heavy particle and the path z{t) of the light one. 

We are now facing some fundamental questions: i) Can the term (17) be 
associated with a local potential V , acting on particle 1? ii) What are the 
Lorentz properties of this potential — does it transform as a scalar or as a 
vector? 

To answer question i) we follow the arguments given in [2]. We must 
therefore return to the exponential function in (3), defining the dynamics of 
the system. The integral in the exponential function in (3) is over the surface 
inside the quark and antiquark trajectories; the characteristic length of these 
trajectories is Tg, a period of quark orbiting. 

Being at some point z{tq) on the trajectory, a quark is influenced by the 
fields and through them by its partner. The radius of nonlocality of the fields 
is given by the correlation length T^, defining behavior of {F{u)F{u')), i.e., 
by the functions D{{u — u')/Tg), Di{{u — u')/Tg). 

Therefore the criterion of local dynamics is Tq ^ Tg [2,14]. 

In the opposite case, Tq <Tg, the quark "feels" all the fields and also the 
motion of the antiquark during all its history. This is the nonlocal dynamics, 
which can be treated by the QCD sum rules [15]. Lattice calculations yield 
Tg ~ 0.2-=-0.3 fm [16], while Tq for both the light and heavy quarks is Tg > 1 
fm. Therefore the actual situation is close to the local dynamics. 

Regarding the first point, we proceed as in [8, 17], forming the minimal 
surface via the connection of 7 = f/s = t/s = t/T, where t = Z4 — y^. 

In this case, combining all exponential functions in (13) and (17), we 
obtain 

T r_2 „ „ 1 , 1 

(18) 



T 


r 1 


B= j dt 



. ^ 


where we have defined, as in [8,18], 


w^{t) = z^{t)p + z^{t){l-p), W^ = —W^, 




r^{t)=z,{t)-z,{t), 2/ii = T/.. 



(19) 



Taking (5) into account, we have 

Wi{t) = Zi{t)l3 + Ri{l- l3), Ri = const, 

w^{t) = z,{t)P+{y^ + t){l-P), (20) 

w^{t) = z,{t)(3 = ri{t)p, W4 = Ziit)P + l-p. 

Following the same procedure as in [8,18], we obtain from B (considered 
as an action) the proper-time Hamiltonian, which depends on the dynamical 

mass /ii: 

The eigenvalue of H{ni) is to be minimized with respect to fii; the final 
value of Hi is determined in this manner. 

Thus we see that the nonperturbative correlator {{F F)) nonpert does indeed 
yield a potential-type term in the proper-time Hamiltonian, provided Tg is 
small [2] and the expansion of the square-root term in f3 is done as in [18] 
(the corresponding error is ~10% ). 

It is more difficult to answer the other question — about the Lorentz 
nature of this potential — since Eq. (21) is written in the cm. system, 
and we do not know how this expression transforms under Lorentz boosts. 
To get a partial understanding of this point, let us compare how the vector 
and scalar potentials enter the Feynman-Schwinger representation for the 
Green's function of a scalar quark [i.e., neglecting a spin term ~ a^^' in (13)]. 

For the scalar case we can write 



G{x,y) ~ / dsDzexp — / Vs{z{T))dr — m^s — / —dr 
while for the vector case we have 







(22) 



G{x,y) ~ / dsDzexp —ra^s — / —dr + i soVf^dZf^ . (23) 

The main difference lies in the fact that trajectories of the Z-type depicted 
in Fig. 2 (virtual pair creation) give different contributions in the scalar and 
vector case. For the scalar case the proper time intervals Ari, At2, and Ar^ 
are positive and add together, thereby suppressing the pair creation, while 
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in the vector case Az^^ , Az^ , and Az^^ strongly cancel each other, making 
pair creation easy. In other words, for the scalar potential all parts of a 
trajectory add arithmetically, while for the vector potential should take a 
vector sum of all intervals along the trajectory. 

Looking at our prototype of the potential, the last term in (18), one can 
see that in this case one actually adds a scalar quantity adtJ~ (we recall that 
dt is not the fourth component of a vector, but rather the proper time cir, 
which is scalar). Thus we seem to be justified in treating the confining force 
as a scalar and not a vector. 

There are additional arguments in favor of this conclusion. 

First of all, in the nonrelativistic derivation of spin-dependent forces in 
nonperturbative QCD [2], the sign and magnitude of the spin-orbit term 
depend on whether the confining potential is chosen as a scalar or a vector 
[4]. In Appendix B we reproduce the spin-dependent terms obtained in [2], 
for the case of two nonrelativistic quarks. One can notice that the asymp- 
totically dominant contribution comes from V/, which contains only a scalar 
contribution and yields the negative coefficient of the spin-orbit term LS, 
characteristic of Thomas precession. 

Results obtained in [2] through the use of the cluster expansion and the 
area law unambigiously predict the spin-orbit force corresponding to the 
Thomas precession term, which was also introduce a previously using the 
string picture [19]. All that corresponds to scalar confinement. 

Second, there are independent arguments in [3] which lead to inequalities 
which are satisfied for the scalar confining potential, and not satisfied for vec- 
tor case. Finally, the vector confining potential is believed to cause the Klein 
paradox [10,11] . In Sec. 4 and 6 we study both scalar and vector confining 
potentials and show analytically and numerically that, indeed, only in the 
first case does one obtain a physically consistent spectrum, corresponding 
to confinement while in the second case one has only quasistationary states. 
This situation is related to the Klein paradox [10,11]. Thus we give strong 
arguments in favor of the scalar confining potential. Unfortunately at the 
moment we are still unable to derive from the Feynman-Schwinger represen- 
tation the Dirac equation with the scalar potential, in the way we have done 
for a perturbative Coulomb interaction. 

Therefore in the next sections we simply postulate the Dirac equation 
with a confining potential of the scalar type (or vector type — to check its 
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inconsistency). 

4 Properties of the spectrum of the Dirac 
equation 

As we showed in the previous sections, the color Coulomb interaction 
between quarks is of a vector nature. We have argued that the confining 
interaction is a scalar. It is instructive, nevertheless, to consider both possi- 
bilities for the interactions, scalar U and vector V, and to study the properties 
of solutions of the Dirac equation in these cases. 

Assuming spherical symmetry, we look for solutions of the Dirac equation 
in the form 

where / + /' = 2j, and the il's are spherical spinors [20]. 
Equations for radial wave functions are 

^ + -G-[E + m + Uir) - Vir)]F = 0, 
dr r 

^--F+[E-m-U(r)-V(r)]G = 0, (25) 

dr r 

where E and m are the energy and mass of the light particle (m = rrii in the 
notation of Sec. 2), and 



(26) 



r -(/ + !), if j = / + l/2, 
1 /, if J = / - 1/2, 

so that \k,\ = j + 1/2 = 1, 2, ... . We consider three choices. 
a) Assume 

U{r) = M\ V{r) = -^, M^ = a; C = ^«s. (27) 

Putting m = and introducing the dimensionless variables x = Mr and 
e = E/M, we arrive at 

G" + -G - I 5 + a; + ^ ) F = 0, (28) 

X \ X J 
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F'--F+(e-x + ^]G = 
X \ X J 

(hereafter the prime denotes the derivative with respect to x). 

Let us find the asymptotics of the radial functions at zero and infinity. 
For a; — >■ 0, inserting the expansions 

G = Ax'' + ..., F = Bx'' + ... 

into (28), we obtain 

iy={K^-CY^ A/B = k-\k^-CY'-'^- (29) 

The wave function is regular at zero for ( <\ k, \= j + 1/2 while for 
C > j + 1/2 there occurs a "collapse to the center" ( well known in quantum 
mechanics ) [21,22]. The difficulty is actually of a formal character and can be 
removed by introducing a cutoff of the Coulomb potential at small distances, 
r < vq — in the same way as in QED with the charge Z > a^^ = 137 (see 
refs. [9-11, 23,24]). 

The behavior of the wave functions at infinity is more complicated. From 
(28) it follows that G, F ~ exp(±a;^/2) as x — > cxo , and the solution with a 
positive exponent is clearly not admissible. 

We accordingly write G, F as 

G ^ A'exp\-{ax^ + bx)]x^ (l-- + % + ..) , 

\ X x^ ) 

F^5'exp[-(ax2 + 6x)]x^ [l- — + ^ + ...J . (30) 

Inserting these expansions into (28) we find, after some elementary but 
cumbersome calculations, 

N = -B', a = 1/2, 6 = 0, P = e^/2, 

c, = {C- 1/2)6, c[ = {C + 1/2)6,.... (31) 

The following coefficients, C2, Cj, ... , can be calculated by means of recurrence 
relationscJ. 



^^Thanks are due to D.Popov for the derivation of the relations and for verification of 
Eqs. (29)-(31). 
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Thus the wave functions fall off at infinity mostly in the same way as 
in the case of the harmonic oscillator in nonrelativistic quantum mechanics. 
The solution satisfying the conditions obtained above, (29) and (30), exists 
only for certain values oi e = SniC, k,), and the wave function is normalizable. 
Therefore, for the Dirac equation only a discrete spectrum exists. 

This conclusion holds for any scalar potential which grows at infinity. 



G{r), F{r) ~ exp 



U{r)dr' 



oo. (32) 



Even for the logarithmic potential the wave functions decrease asymptot- 
ically faster than exponentially: 

G,F(_r~*oo) ~ exp(— (?rlnr) for U{r) = glnr. 

b) When both interactions are of a vector type, i.e., 

U{r) = 0, V{r) = -Cr'^ + M^r, (33) 

the system of radial equations is 

G' + -G-Ie + --x]f = 0, 

X \ X J 

F'--F+(e + ^-x]G = 0. (34) 

X \ X J 

In this case, at a; —> we have 

while the parameters in Eq. (30) assume the values 

A' = iB', a = -i/2, b = ie, f3 = -<• 

Here we chose a solution satisfying the Sommerfeld radiation condition, 
i.e., having an outgoing wave at r — i> cxd. Thus for interaction (33) the discrete 
spectrum is absent (see also curve 2 in Fig. 3), and all solutions of the Dirac 
equation, if any, are quasistationary. In the special case of a square well one 
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can obtain an exact analytic equation for the spectrum [25]. The width of 
these quasistationary states determines the probabihty of spontaneous pair 
creation in potential (33). This phenomenon is closely connected with the 
Klein paradox [10,11]. Physically it means that the problem is now of the 
many-particle type, but actually the Dirac equation is still applicable. 

c) For completeness, consider also the case of purely scalar interactions, 

f/ = -(r"i + MV, V = 0. (35) 

In this case we have at r ^ 0, instead of (29), 

u={n' + CY'\ A/B =[n- {k' + CY"]K. (36) 

Thus the "collapse to the center" does not occurs for arbitrary large 
values of the Coulomb parameter C,. At infinity we obtain a behavior which 
is analogous to (30): 

G,F~exp(-MV/2); 

therefore the spectrum is discrete. 

Hence the character of the energy spectrum depends crucially on the 
assumption about the nature of the confining interaction M'^r (scalar or 
vector). 

The reason can be understood qualitatively by using the method of an 
effective potential [24,26]. 

System (25) corresponds in the quasiclassical approximation (neglecting 
spin-orbit and spin-spin forces) to the following relation between energy and 
momentum 

p^ir) = [E- V{r)f - [m + U{r)f + /^V'^ = 2m{E - Uejf), (37) 

where E and f/e// are the effective energy and potential in the nonrelativistic 
Schrodinger equation: 

E = —{E^-m'), U,ff = —{U''-V^) + U + ^. (38) 

2m^ '' ^^ 2m^ ' 2mr^ ^ ' 

In the nonrelativistic limit {E ^ m and f/, | V |<^ m) we have 

Ueff^U + V + {l + l/2f/2mr^. 
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Therefore, all three cases, a)-c) correspond to the tunnel potential 

Ueffir) = -C/r + ar, 

which is often used in QCD. 

However, for large U and V, the situation is drastically different for our 
cases a)-c). This can be understood from Fig. 3, where the qualitative 
behavior of Ueff{r) is shown. Using the effective-potential method, one can 
also perform a quantitative study of the problem (for details see Appendix 
A). Applying the WKB method and taking into account Eq. (37), one can 
easily determine the asymptotics of the wave functions at infinity. 

For example, in case a) we get 

-jWxU \p\ = U + m-^+^— + .... (39) 

If U{x) = gx"", a > 0, and V{x) = ~C/x, then at a; — ;► oo we find 

_ j exp{-[gx^/2 + mx-{2g)-^£^\nx]}, a = l, , . 

^ \ exp{-[g{a + iy^x''+^+mx]}, a>l. ^ ' 

The first formula (with g = 1 and m = 0) explains the structure of 
asymptotic expansion (30). 

5 Symmetry properties of the spectrum of 
the Dirac equation 

Here we discuss symmetries of solutions of the Dirac equation for particles 
of zero mass. 

It is clear that when both m and U are zero the Dirac equation (and the 
corresponding term in the Lagrangian) is chirally symmetric, i.e., does not 
change under a global transformation: 

^ -^ exp(m75)^, ^ -^ 'I'exp(ia75) . (41) 

From the point of view of the spectrum the chiral symetry means that all 
states are parity degenerate; i.e., the masses of the states 0"*" and 0~ (or l"*" 
and 1~) are the same. 
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It is easy to demonstrate that system (25) with m = is invariant under 
even more general transformations when U ^ 0: 

E^E, K-^-K, U -^-U, V ^V, 

G{r) -^ -F{r), F{r) -^ G{r). (42) 

It follows from (42) that for f/ = the spectrum is degenerate with respect 
to the sign of the Dirac quantum number k\ i.e., it depends on only the total 
momentum j, not on the way / and s are coupled (chiral degeneracy). 

Another symmetry of the Dirac equation is 

E^-E, K^-K, U ^U, V ^ -V, G{r) ^ F(r), (43) 

which in contrast to (42) connects states with positive and negative energy. 
In particular, for V = (scalar interaction) there is a doubling of states with 
a given \ En \,En = ^ \ En \, and we can always consider En > 0. In the 
specific case E = (zero modes) there is, at first sight, a chiral degeneracy 
of states. However, it is easy to show that the degenerate states with E = 
do not exist at all. Multiplying the first equation in (25) by G, multiplying 
the second by F, and integrating the difference from r = to r = oo, we 
obtain an identity: 

oo oo 

f {G^ + F^)— = 2k'^ f {E - V)GFdr (44) 



(we have used in the derivation the circumstanse that G and F vanish at 
both r = and r = cxd). From Eq. (44) one can see that for the purely scalar 
interaction there are no solutions with zero energy. 

Note that the symmetry transformation, Eq. (42), can be also presented 
in an operator form. For a zero-mass particle the Hamiltonian H and the 
Dirac operator K are 

H = ap + [3U{r) + V{r), K = -(3{(t\ + 1) (45) 

(the quantum number n is the eigenvalue of the operator K). It can be easily 
seen that in this case we have 

75^^75 = ap - Wr) + V(r) , 75^75 = "^^, (46) 

which are the same as Eqs. (42). 

16 



6 Results of numarical calculations 

Here we report results of our numerical calculations. In Fig.4 is shown the 
dependence of the eigenvalues £„« = Enn/M on the ratio (/ \ n \ioY potential 
(27) ( the solid lines correspond to the lowest states for a given value of k, 
n = 1; the dashed lines are for the first excited states, n = 2). The energy 
eigenvalues decrease monotonically with growing Coulomb parameter (, and 
for (^ — >| K I they have a square- root singularity. The latter is characteristic 
of potentials with a Coulombic behavior at r — > and is connected to the 
phenomenon of the collapse to the center [10,11]. 

In Fig. 4 one can see that the levels with k > lie much higher than 
those with k < (for a given n). The physical meaning of this effect becomes 
clear when one recalls that the centrifugal barrier is proportional to k,{k,+ 1) 
and is absent for k = — 1 (for example), in contrast to the case oi k = +1. 

Note that the energies of the lowest (n = 1) states with k < reach the 
value e = at the maximal possible value of the Coulomb coupling constant 
( = —K (Fig. 5). All the other states also have square- root singularities 
at (^ = |/t;|, but their energies are positive (Fig. 6). The numerical values of 
^riHiiC) ill the two extreme cases C = and C = 1^1 are given in Table I, where 
the dependence of Suk. on the quark mass is also illustrated. Figure 7 shows 
the energies of several lowest levels within the interval 0.3 < C < 0-8; which 
is of practical interest for the vh system ell. 

We note that the singularity at C = |fi;| can be removed when one intro- 
duces a cutoff of the Coulomb potential V{r) = — C/^ at small distanced^: 

^W = lr^;/\.r / . 'J n^'°' (47) 

1 (-C/^o)/(r/ro), if < r < ro, ^ ^ 

where ro is the cutoff radius, and /(O) < oo. 



^■'As can be seen from Figs. 5 and 6, the dependence of the energy e„K on the hght-quark 
mass is practically linear. This is in agreement with ref. [27], where the energy spectrum 
of Eq. (25) was calculated at C = 0.6 with hyperfine splitting taken into account, which is 
necessary for a detailed comparison with the experimental spectrum (the procedure was 
similar to the calculation of the hyperfine structure for relativistic Coulomb problem). 
The value of hyperfine splitting obtained in ref. [27] differs from experimental data for 
charmonium and bottomonium by over 30%. 

■^^In the case of the lowest levels with k < 0, for which £(k) -^ in the limit C ^| k | 
(see Fig. 4), this singularity can be derived analytically — see Eq. (50) below. 
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In the case tq -C M~^ the influence of the cutoff on the energy levels is 
important only in the immediate vicinity of the point ( =\ k \. Here the 
level sinks to the boundary e = (it corresponds to a boundary e = — m for 
m 7^ 0) at some value ( = (nn >| '^ I5 which depends both on Tq and on the 
concrete form of the cutoff function /(r/ro) in Eq. (47). 

With a further increase of ( the level goes on down to the region of 
negative values of energy, but the spectrum stays discrete (this is in contrast 
to the Coulombic problem with a vector potential [10,11]). Thus pair creation 
does not take place. 

7 Exact solution for E = 

In QED the solutions of the Dirac equation simplify considerably for 
E = ±m, which corresponds in our case (m = 0) to ii^ = 0. Eqs. (28) in 
terms of linear combinations u = G + F, U = G — F assume the form {e = 0) 

u' — xu -\ V = 0, v' + XV -\ u = 0. (48) 

X X 

At C = —fi the first equation can be solved explicitly: 

u = Cexp(a;^/2), 

which yields C = and hence G = —F = v/2. The normalized wave 
functions areell 

G{x) = -F{x) = 7r-^/^exp(-a;V2). (49) 

Using the perturbation theory in the parameter (/t^ — C^)^^^ ^ 1; one can 
determine the behavior of the energy near ( = —k: 

.(C) = vr-ni-CV'^Y/' + - (50) 

(for details see Appendix A). Computations confirm this asymptotic expan- 
sion; see Fig. 5. 



■*^This result can be easily generalized to the case to ^ and arbitrary scalar potential 
U{r). 



Note that the solution of the type in (49) exists only for the states with 
K < 0. In the case C = '^^ Eqs. (48) have a nonzero solution 

G = F = const • exp(a;^/2) 

which is not admissible because of the exponential growth at infinity. 

These results can be generalized to arbitrary values of ( and n. Solving 
(48) for one of functions u, v, we come to the equation 

,2 /■2' 

w = 0, (51) 



w H — w — 

X 






x^ 



where the upper (lower) sign corresponds to the function u (f ). We consider 
the case C >| ^^ |, when the cutoff of the Coulomb part is essential. A solution 
of these equations for r > tq, decreasing at infinity, can be expressed in terms 
of Whittaker functions, 

U = CiX-'W-i/2,igix^), V = C2X~'Wy2,g{x^), (52) 

where g = (C^ — k'^Y^'^ > 0, and Ci, C2 are constants. 

Insertion of (52) into the first of Eqs. (48) yields the connection between 

Cl,C2'- 

Ci=2{k + C)C2. (53) 

In the internal region, < r < tq, the Dirac equation should be solved 
with cutoff potential (47), where one can neglect the linear potential because 
of the relation tq <^ 1/M. For the simplest case, f{x) = 1 (a square-well 
cutoff), the solution can be found analytically and can be expressed in terms 
of Bessel functions with half-integer index. 

The energy spectrum is found from the joining of the internal and external 
solutions at r = tq. In the case of the Dirac equation one can join the ratio 
F/G instead of the logarithmic derivative. 

As a result we have 

2{n + OW_y2,^g{xl)/W^/2,^g{xl) = ^, (54) 

where ^ = u{xo)/v{xq) is determined from the interval solution and depends 
on K, ( and the cutoff model. For the states with k, = =1=1 in the square-well 
cutoff case, we have 

n - 1 i-C(i + ctgC) .... 

^^^^-'~ e-^"i + C(i-ctgC)- ^^^^ 
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Equation (54) can be solved numerically; it determines the dependence of 
the Coulomb constant (nn, which corresponds to the descent of the n-th level 
to the boundary e = 0, on the cutoff radius tq [for a certain choice of the 
cutoff function /(r/rg)]. 

Finally, we note that a solution (49) can be generalized to the case of 
nonzero quark mass and an arbitrary scalar potential U{r). There is an 
exact solution of the Dirac equation with ( = —k and i? = 0, 



G'o(r) = -Fo(r) = iVexp 



m{r')dr' 







(56) 



where m{r) = m + U{r) is a variable quark mass and A^ is the normalization 
constant, [see Eq. (A. 8)]. Since this wave function has no nodes, this solution 
corresponds to the lowest (n = 1) level with the fixed k. The energy of the 
state, as well as other physical quantities, has a square-root singularity at 
(^ -^ —K, characteristic of the relativistic Coulomb problem: 

E{0 = c^P + O{P^), p=l-c,P+..., (57) 

where P = (1 — C^//t^)^/^, and p is the relative weight of the lower component 
of the Dirac bispinor, 

oo oo 

p = I F\r)dr / f G\r)dr. (58) 







It can be shown that 



ci=N^ , C2 = 4N^{r) (59) 

(see Appendix A), where (r) is the mean radius of the bound state (56). The 
effective-potential method given in Appendix A is most useful in deriving 
these formulas, as well as for a quantitative analysis of the wave functions 
near ( = ±k. 

As is seen from Eq. (57), at C ~ \i^\ the parameter p is close to one 
(see also Table II). This means that the motion of a light quark is definitely 
relativistic. The coefficients Ci and C2 for a particular case oi m{r) = m + crr 
are presented in Fig. 8. Finally, let us note that square-root singularities 
of the type of (57) are directly connected with the behavior of the Coulomb 
potential at r ^ and disappear when it is regularized at small distances. 
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8 Discussion and conclusion 

In conclusion we discuss the structure of the resulting spectrum, briefly 
compare it to experiment, and summarize our resultsE^. 

We start with the nonrelativistic spectrum and use the formulas given 
in Appendix B, which yield the spectrum shown in Fig. 9, a. Here the 
splitting between the 5" = 1 and S" = (L = 0) levels is due to the hyperfine 
interaction and, as is seen in Appendix B, is proportional to {mim2)~^ . When 
both masses are large {rrii ^ M = cj^^"^)., all splittings are small, including 
the spin-spin [0((mir?7,2)~^)] and spin-orbit splittings, which contains terms 
0{mi'^), 0(1712'^) and 0{{mim2)~^), which one can denote as {LS)i, {LS)2 
and (LS) 12 — see Appendix B. 

In the case m2 3> mi ~ M (Fig. 9,6), the intervals of {LS)i become 
the largest, and the spectrum transforms in such a way that the coupling of 
the spin of particle 2 becomes very weak. In the limit m2 ^ 00 it finally 
decouples, leading to the Isgur-Wise symmetry [5]. Namely, the states can 
be classified by the total momentum of the light particle j, while the states 
of total momentum of the system J = j + S2 are degenerate with respect to 
the direction of 82- 

This is what one observes in the Dirac spectrum (see Fig. 4 of our nu- 
merical calculations and Fig. 9, b, which give a schematic description of the 
Dirac spectrum). Here the degeneracy (Isgur-Wise symmetry) is complete. 

This should be compared to the experimental picture in Fig. 9, c, where 
for the charmed mesons the order of the levels is the same as in Fig. 9, b , 
but the splittings are still large. For 5-mesons the splitting of lowest levels 
is smaller by a factor of 3 (52 MeV), as it should be, since the mass of the 
c quark is ~ 3 times smaller than that of b quark. We also note that the 
experimental splitting AE ~ 450 MeV between j = 1/2 and j = 3/2 states is 
reproducible in our results when one takes ^/a ~ 0.5 GeV and ( ~ 0.6 4- 0.8. 

Thus we can conclude that the Dirac equation yields a reasonable quali- 
tative description of the actual spectrum. We hope to discuss this point in 
more detail in future publications, where we will also give predictions for the 
Ds-, B-, and i?s-mesons. 

Summarizing our results, we have derived the Dirac equation from the 
Feynman-Schwinger representation of the quark-antiquark Green's function 



^^The main results of this paper were presented in ref. [8]. 
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in case of the color Coulomb interaction. 

Assuming that the confining interaction can be introduced into the Dirac 
equation in the sameway as to the Coulomb interaction, we have clarified the 
Lorentz structure of the confining interaction connecting a light quark to a 
heavy antiquark. 

The analysis of the solutions of the Dirac equation shows that a potential 
growing at infinity yields confinement only if it is a scalar, not the time 
component of a 4-vector. In this aspect there is an essential difference from 
the nonrelativistic case. 

We have also studied the dependence of the charge ( in the critical region, 
C ~ |k|, and found its dependence on the cutoff of the Coulomb potential. We 
have computed energy eigenvalues for several lowest levels and have compared 
them with the nonrelativistic spectrum and experimental results. 

After this paper had been finished and submitted for publication, we 
learned of several papers [28-30] in which analogous problems are discussecCll. 

In the paper by Ono [29] , a relativistic generalization of the Richardson 
potential was considered. The Coulomb part of the potential was considered 
as a 4-vector component, and its confining part as a relativistic scalar. The 
paper by Dremin and Leonidov [30] discusses the properties of solutions of 
the Dirac equation with scalar and vector coupling. Using the squared Dirac 
equation, the authors arrived at the conclusion that it is only in the case 
of a scalar potential growing at infinity that bound states can exist in the 
system. That conclusion is in full agreement with the results of our Sec. 4, 
where we study in greater detail the asymptotics of the Dirac wave functions 
at infinity and give a physical explanation of the above results using the 
effective potential method. 

APPENDIX A 
Method of effective potential 

The system of Eqs. (25) can be reduced to an equivalent Schrodinger 
equation with a potential depending on energy and having, in general, a 
rather complicated form O. There is a considerable simplification if (^ = 
\k\ = 1,2 This case we shall discuss further. 



^^The authors are indebted to A.E.Kudryavtsev who pointed out to us the above men- 
tioned papers. 

^^See refs. [10,24,26] on physical phenomena near the boundary of the lower continuum 
E — —mc^ in QED [10], where the above method was successfully applied. 
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First let ( = —k. Proceeding in Eqs. (27) to the linear combinations 
u = G + F^ V = G + F, we obtain equations from which the function v is 
easily excluded. Then we get 



u" +[e^ -l-{x + fiy + 2Cex-^]u = 0, u{0) = 0, 

l + ^]u\x)dx = i (A.l) 

^ ^^^ 

i.e., the Schrodinger equation with the effective energy E and potential U 

E = 1(5^-1), t/=i(x + /if-eCx^^ (A.2) 

Note that the boundary condition at zero for u{x) follows from the fact that 
we have A/B = —1 in Eq. (29), and the normalization condition in (A.l) 
corresponds to the usual condition 



f{G^ + F^)dr = 1 



for the discrete spectrum. 

In particular, for e = we have E = —1/2, which cannot be an eigenvalue 
for the harmonic oscillator. Therefore, we have u{x) = and v'+{x+iJ,)v = 0, 
from which the exact solution (49) follows. 

For the case ( = k, the function v, but not u, satisfies an equation of the 
like (A.l), and we have E = {e'^ + l)/2. Here we have E = 1/2 with e = 0; 
therefore we have v{x) =0 (with the boundary condition at zero taken into 
account), while all nonzero solutions u{x) increase exponentially at infinity. 
Thus, the Dirac equation has no physically acceptable solutions with zero 
energy ii ( = k, > 0. 

The following results are obtained from (A.2). Let En = En{<y,fi) be the 
energy spectrum of the Schrodinger equation with the potential 

V{r) = -{r + fif-ar-\ (A.3) 

with < r < oo and / = [in particular, i5„(0,0) = 2(n — 1/4), n = 
1,2,3,...] . Then the eigenvalues e^^ of the problem considered may be 
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obtained from the transcendental equation: 

E„(|/t|e,/i) = -(5^±sgn k) 



(A.4) 



(the ± sign coincides with the sign of the quantum number k,). Due to the 

relations 

dEr,. /1\ dE„. , . , 



da \r/ ' dfi 

sometimes called the Hellmann-Feynman theorem, the energies En{<y,fi) de- 
crease monotonously as the Coulomb parameter a goes up and increase with 
the growing /i. Taking this into account, we get 



4"^ = o<4^)<5«<4^)<... 



(A.5) 



(with fixed |/t|). The results of the numerical solution are given in Table I. 

The results are easily generalized to the case of an arbitrary scalar po- 
tential U{r). Instead of (A.l) we get 






E' 



[m + U{r)f 



dU 2\k\E' 



^sgn.j- + 



x = o 



(A.6) 



where x(0) = 0, and the normalization condition does not depend explicitly 
on the form of U{r): 

|/t| \ o, . . 1 



^ + '^Jx\r)dT ^ 



(A.7) 



[x = '"(^) for the states with (^ = —k, and x = '^(^) for those with ( = k]. 

At ( = —K, there is an exact solution of the Dirac equation which corre- 
sponds to zero energy [and does not depend explicitly on k; see Eqs. (54)]; 
here Xo(^) = 0, 



G'o(r) = -Fo(r) = A^exp 



m{r')dr' 



N = (2 f drexp -2 f m{r')dr' 



-1/2 



(Ai 
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The behavior of the energy of the state at ( -^ —k, can be also determined. 
Assuming 



u = I3ui{r) + ... , V = vo + (3vi{r) + ... , 
and substituting these expansions into (25), we obtain 

dui ^ dvi 2\k\ 

— mui + Cifo = 0, — — h mvi = Ui. 

dr dr r 



(A.9) 



(A.IO) 



The solution of the first equation which decreases at infinity has the form 



■ui(r) = 2ciA^exp 



m{r')dr' 



(ia;exp 



2 / m{x')dx' 



(A.ll) 



while the function fi(r) is calculated by a quadrature; at r = we have 
Mi(0) = ci/iV, t;o(0) = 2N. Let 

then it follows from Eq. (28) that 

^^' A-B (l + /?)i/2 + (1-/5)1/2 2'^'^^l' > 

whence Mi/vo|r=o = 1/2, and finally we get Eq. (59) for Ci. 

Thus for the states whose energy vanishes at C = ~i'^) the coefficient of 
the square-root singularity in E{C) can be found explicitly for an arbitrary 
scalar potential U{r). Other physical quantities, such as p, have the same 
singularities at C = ~'t- 

Using Eqs. (25) and (A. 6), we find for the parameter p introduced in 
(58) that 

oo oo 

= 2 / uvdr / [u^ + v'^)dr = 





1 + s' 



oo oo / I I \ ,1 °^ 

m(r)x^dr / i E -\ X^dr = — / m{r)x^dr 

^ ^ 



(A.12) 
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oo 


r 


C2 


= 8cl / exp 


-2 fm{r')dr' 




J 




J 





[here C = \i^\, and in the latter formula x(^) is normalized according to Eq. 
(A. 7)]. Hence, expansion (57) for p follows, where 

oo 

C2 = 2 / uiVodr. 



Substituting in the above expressions for ui,vo, and changing the order of 
integration, we obtain 



rdr, (A. 13) 



which coincides with Eq. (59). 

Let us consider a few simple dependences m{r). 

1) A constant mass m{r) = m corresponds to QED. At C = ^'tj we have 

Go = -Fo = m^/^e-™", Ci = m, Cs = 2 , (A.14) 

which corresponds to the Sommerfeld formula for the relativistic Coulomb 
problem, which has the following form for the k, = —n states (j = / + 1/2 = 
n-1/2): 

/ CV^^ 1-0 

En.iC) = m ^l - ^j ^mP, P=Y^- (A.15) 

The levels with radial quantum numbers n,. > also have the Coulomb 
singularity at C, — * |/t|, but the corresponding energies are positive. For 
example, for the nsi/2 states {k = —l,nr = n — 1): 



/ ji — \ R \ /2n 

(A.16) 
where (3 = {I - Cf/^. C = Za^l.anAN = (n^ - 2n + 2)^/2. 

At n ^ 1, we find p„ = l/4n^ — > 0, and the bound states become 
nonrelativistic. 

2) At m{r) = m + ar, we find 

Go = -Fo = Nexp{-{ar^/2 + mr)}, (A.17) 
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N = (a/7r)i/^[exp(/i2) . erf c(/i)]"^/2^ fi = mja^l'', 

ci = iV^, C2 = 4a-^ci(ci - m) (A.18) 

(note that the dimensionless coefficient C2 depends only on yu). In the extreme 
case at /i ^ 1 we obtain 

'' = G)'" (' + ^75 + ■■■) ■ -= = 7 (1 + 1^7?" + ■■■) • (A-19) 

while at /i ^ 1 we obtain 

ci =m(l + -/i-2 + ...), c2 = 2-/i-2 + ... (A.20) 



Numerical results for ci, C2 are given in Fig. 9. 



APPENDIX B 



The nonrelativistic spectrum of two particles with all spin interactions 
taken into account can be obtained in the perturbative way {1/m expansion): 

/ I /cTiL cr2L\ ( 1 de 1 dVi\ 
m = nil + m2 + E{n) + (n 



nil "^2 / \4^ dr 2r dr 
(cTi + cr2)L 1 dV2 (Tia2 



2niini2 r dr 12niini2 

+ —- (ScTin-CTan- cricr2)V'3(r) I n). (B.l) 

llniini2 ' 

The contributions of the scalar interaction U{r) = ar and the vector in- 
teraction V{r) = —4as/3r to the spin splittings are 

1 de a 4 as 1 , u 1 , Aa^ 

r dr r 6 r"^ r r r 6 r'^ 



V.^?^Si^>(r). V3 = % (B.2) 

6 r-^ 

The main difference appears in VI, where only the scalar interaction en- 
ters. Note that the overall signs of the spin-orbit terms due to the scalar 
(ar) and vector (— 4Q;s/3r) potentials are different. 
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TABLE I. The eigenvalues Euk of the Dirac equation with C = and ( 



\K\ 





C = 


= 


c = 


\k,\ 


n 


K 


/i = 


/i = 0.3 


fi = 


/i = 0.3 


1 


-1 


1.61944 


1.84441 


.00000 


.00000 


2 


-1 


2.60263 


2.80689 


1.06901 


1.17491 


3 


-1 


3.29118 


3.49080 


1.96846 


2.11712 


4 


-1 


3.85541 


4.05300 


2.67856 


2.84277 


1 


-2 


2.14652 


2.36761 


.00000 


.00000 


2 


-2 


2.95197 


3.15853 


.64443 


.70949 


3 


-2 


3.57353 


3.77508 


1.39768 


1.51077 


4 


-2 


4.09947 


4.29854 


2.08531 


2.22342 


1 


-3 


2.56927 


2.78850 


.00000 


.00000 


2 


-3 


3.26852 


3.47647 


.45001 


.49567 


1 


-4 


2.93218 


3.15029 


.00000 


.00000 


1 


1 


2.29403 


2.49206 


.64015 


.79011 


2 


1 


3.03103 


3.22747 


1.62588 


1.79736 


3 


1 


3.62598 


3.82161 


2.39019 


2.57021 


1 


2 


2.70440 


2.90645 


.36916 


.46149 


2 


2 


3.35376 


3.55217 


1.12003 


1.25059 


1 


3 


3.05967 


3.26183 


.25449 


.31931 


1 


4 


3.40866 


3.60000 


.19328 


.24284 
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TABLE II. Values of the parameter p in Eq. (5^ 





C 


= 


c = 


K 


n 


K 


/i = 


;u = 0.3 


p = {) 


;u = 0.3 


1 


-1 


.15534 


.13156 


1.00000 


1.00000 


2 


-1 


.20341 


.17727 


.50983 


.44930 


3 


-1 


.21280 


.18976 


.35668 


.31909 


4 


-1 


.21640 


.19580 


.30697 


.27815 


1 


-2 


.16165 


.14188 


1.00000 


1.00000 


2 


-2 


.19556 


.17401 


.66958 


.61953 


3 


-2 


.20662 


.18651 


.47148 


.43447 


4 


-2 


.21178 


.19318 


.38478 


.35496 


1 


-3 


.16444 


.14718 


1.00000 


1.00000 


2 


-4 


.16601 


.15049 


1.00000 


1.00000 


1 


1 


.22266 


.18839 


.34861 


.31937 


1 


2 


.20986 


.18360 


.54374 


.51578 


1 


3 


.20203 


.18024 


.65694 


.63317 


1 


4 


.20532 


.17791 


.72679 


.70663 
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Figure captions. 

Fig. 1. Quark-ant iquark Green's function describing transition from the 
initial state r{y,y) to the final state r{x,x). Wavy lines refer to the quark 
(antiquark) Green's function Ai{x,y) {A2{x,y)) 

Fig. 2. Trajectory corresponding to virtual pair creation, a — Scalar 
case; b — vector case 

Fig. 3. Nonrelativistic potential Vj^/ji{r) (a) and effective potential Ueff{r) 
(b) for various Lorentz structures of the interaction. Curves 1-3 correspond 
to the three choices in Sec. 4 

Fig. 4. Energy spectrum for the Dirac equation with potential (27). The 
solid curves correspond to the n = 1 states, the dashed ones to n = 2. The 
curves are labeled with the values of k, 

Fig. 5. The eigenvalues ei^i^O, n = 1, of the Dirac equation with po- 
tential (27) in the vicinity of ( = —k. The solid curves correspond to the 
K = —1 states, the dotted ones to k = —2. The curves are labeled with the 
values of /i = m/M 

Fig. 6. The same as in the preceding figure, for excited states {k = 
— 1, n = 2 and k = l,n = 2) 

Fig. 7. The energies of the lowest levels, Enn, versus the Coulomb param- 
eter ( (for a quark mass m = 0) 

Fig. 8. Dependence of ci and C2 [Eq. (59)] on the quark mass 

Fig. 9. Energy spectrum for the heavy quark-light antiquark system 
(qualitative), a — mi,m2 ^ M = a^l"^; b — mi ^ 0,m2 — * oo ; c — 
experiment. The energy scale in part a is much larger then in parts b and c 
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